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The  contents  of  this  memo  represent  the  written  version  of 
a  presentation  at  the  Ocean  Se ismic- Acoustic  Conference  held  at 
LaSpezia,  Italy  June  10-14  1985. 


ABSTRACT 


Long-range  sound  propagation  in  the  Arctic  Ocean  is 
characterized  by  a  refractive  surface  sound  channel  with  a 
rough  water-ice  interface.  Recent  experimental  measurements 
show  that,  for  frequencies  below  100  Hz,  attenuation  exceeds 
sea  water  absorption  by  two  orders  of  magnitude.  The  most 
likely  mechanism  is  scattering  from  the  rough  ice  canopy. 
Theoretical  estimates  of  the  scattering  loss,  obtained  using 
the  method  of  small  perturbation  and  statistical  measures  of 
the  under  ice  roughness  obtained  from  experimental  data,  were 
examined  for  models  of  the  ice  canopy  having  varying  degrees  of 
realism.  All  theoretical  estimates  for  scattering  loss, 
irrespective  of  the  particular  model  for  the  ice  canopy,  were 
substantially  lower  than  the  measured  values  of  scattering 
loss.  The  physics  of  the  loss  mechanism  is  evidently  not  well 
understood  and  evidently  additional  experimental  and 
theoretical  investigations  are  required. 
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LOW  FREQUENCY  ATTENUATION  IN  THE  ARCTIC  OCEAN 


The  deep  water  abyssal  plain  regions  of  the  Arctic 
ocean  provide  almost  text  book  examples  of  Acoustic 
Propagation  in  a  range  independent  ocean  at  low  fre¬ 
quencies  if  one  could  ignore  the  influence  of  the  rough 
canopy.  This  can  be  seen  by  examining  the  empirically 
derived  expression  and  plots  of  propagation  loss  versus 
range  shown  in  Fig.  1. 

Shown  at  the  top  is  an  expression  for  propagation 
loss  empirically  derived  from  data  in  an  abyssal  plain 
of  the  Arctic  ocean.  The  first  two  terms  (A  +  101og(R)) 
represent  the  expected  result  for  a  range  independent 
ocean.  The  last  term  («R)  is  similar  in  form  to  the 
term  representing  attenuation  loss  found  for  ice  free 
waters.  We  note  that  for  the  Arctic,  however,  this  term 
is  predominantly  determined  by  scattering  loss  from  the 
rough  ice  canopy,  and  that  it  is  two  orders  of  magnitude 
larger  than  attenuation  loss  in  ice  free  waters. 

The  scattering  loss  was  found  to  have  an  f 3/2 
dependency  which  results  in  a  rather  significant 
increase  in  propagation  loss  at  low  frequencies  and  long 
ranges  as  can  be  seen  from  an  examination  of  the  curves. 

This  paper  summarizes  our  progress  in  an  ongoing 
effort  to  gain  insight  into  this  scattering  problem. 

Measurements  of  the  acoustic  field  and  the  under¬ 
side  ice  profile  were  obtained  concurrently.  A  typical 
section  of  the  ice  profile  is  depicted  in  the  upper  left 
hand  corner  of  Fig.  2.  We  have  tipped  the  ice  upside 
down  so  that  the  air-ice  interface  is  represented  by  the 
x-axis.  This  surface  is  also  rough  but  appears  as  a 
smooth  line  owing  to  the  lack  of  concurrent  topside 
roughness  measurements.  The  wavy  line  above  the  x-axis 
is  the  water-ice  interface.  We  call  your  attention  to 
the  fact  that  this  is  a  highly  compressed  view  since  the 
horizontal  scale  is  in  km  and  the  vertical  scale  is  in 
m.  Although  the  deepest  ice  keel  shown  here  is  only 
about  10m,  ice  keels  on  the  order  of  30m  are  not 
uncommon. 
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From  the  distribution  of  ice  keels  shown  in  the 
upper  right  hand  corner  we  find  that  the  mean  draft  is 
4m  and  the  standard  deviation  of  ice  keels  taken  about 
the  mean  is  2m. 

The  normalized  auto  correlation  function  is 
provided  in  the  lower  right  hand  corner.  Using  the 
definition  that  the  correlation  length,  1,  is  that  value 
of  for  which  the  correlation  function  has  the  value  of 
.6  (or  as  it  turns  out  for  which  it  decays  to  one  over 
e)  we  find  that  1  equals  44m. 

The  resulting  one  dimensional  power  spectral 
density  is  shown  in  the  lower  left  hand  corner  along 
with  a  smooth  curve  obtained  from  an  analytic  function 
to  be  described  shortly.  The  spectrum  is  found  to  have 
F  to  the  minus  three  frequency  dependence. 

The  values  for  the  parameters  shown,  namely  a  mean 
draft  of  4m,  -=2m  and  a  correlation  length  of  44m  are 
the  average  values  for  these  parameters  resulting  from 
concurrent  measurements  of  the  ice  profile  and  acoustic 
field  obtained  over  a  large  portion  of  the  abyssal  plain 
where  the  experiment  was  conducted. 

The  availability  of  concurrent  measurements,  which 
is  a  rather  rare  occurence  in  acoustic  scattering,  has 
caused  us  a  fair  amount  of  frustration  in  our  attempt  to 
understand  the  scattering  process.  This  frustration 
stems  from  our  inability  to  obtain  the  correct  answer 
when  using  the  average  values  obtained  from  the 
concurrent  measurements  and  a  theory  which  we  believe  is 
appropriate.  On  the  other  hand,  if  the  concurrent 
measurements  are  ignored,  which  is  equivalent  to  using 
improper  values  for  free  ice  statistical  parameters,  and 
the  same  theoretical  development  is  employed,  it  is 
exceedingly  easy  to  obtain  an  answer  which  appears  to  be 
correct.  The  remainder  of  the  paper  provides  an 
elaboration  of  these  comments.  Our  progress  to  date  can 
thus  be  summarized  as  follows:  we  have  not  found  a 
satisfactory  solution,  but  we  can  indicate  what  the 
answer  cannot  be. 

Candidate  mathematical  models  of  the  ice  canopy  are 
arranged  from  top  to  bottom  in  decreasing  order  of 
realism  in  Fig.  3.  Model  I  being  a  full  elastic  layer 
is  perhaps  the  most  realistic  description  but  it  also 
has  the  greatest  number  of  free  parameters.  Only  the 
density,  the  compressional  sound  speed,  (Cj'p,  the 

standard  deviation  of  the  bottom  roughness,  'g  ,  'and  the 
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correlation  length  of  the  bottom  roughness,  lg  >  can  ^>e 
considered  approximately  known. 

Model  II  ignores  the  elastic  properties  of  ice  but 
still  has  three  free  parameters  namely,  <7T,  3^  ,  and 
which  represents  absorption  loss  within  the  ice. 

Finally  Model  III  essentially  ignores  everything 
about  the  ice  except  for  its  roughness  but  has  the 
attractive  feature  that  no  free  parameters  exist  if  the 
scattering  surface  is  taken  as  the  bottom  of  the  ice  and 
measured  values  are  used  for aB  and  1  .  We  thus  start 
with  Model  III  and  work  our  way  up  in  sophistication. 

The  application  of  the  method  of  small 
perturbations  (MSP)  to  underwater  acoustic  scattering 
problems  has  received  considerable  attention  (1-10)  in 
the  past  and  since  we  believe  it  is  applicable  to  this 
problem  it  is  utilized. 

The  general  procedure  outlined  in  Fig.  4  is  to 
utilize  the  method  of  small  perturbations  (MSP)  to 
obtain  an  expression  for  the  coherent  scattered  j^lane 
wave  reflection  coefficient  which  is  denoted  by  R. 

We  make  the  following  observation  regarding 
perturbation  theory.  First  we  believe  it  is  applicable 
at  the  frequencies  of  interest  and  second  that 
experimental  results  support  neglecting  the  incoherent 
component  of  the  scattered  field  for  this  problem. 

We  next  replace  the  unperturbed  reflection 
coefficient  denoted^Jsy  R,  which  appears  in  the  Fast 
Field  Program  with  R,  insert  a  measured  Arctic  sound 
speed  profile  and  compare  FFP  predictions  to  the 
empirical  curves  shown  previously. 

In  utilizing  MSP  we  make  the  assumptions  that  the 
statistics  of  the  rough  scattering  surface  are 
stationary  and  isotropic.  Our  measurements  of  the  ice 
statistics  along  a  line  indicate  that  the  statistics  are 
not  stationary  and  if  the  topside  and  bottoms ide 
roughness  characteristics  are  even  loosely  correlated 
then  observations  of  the  two  dimensional  topside  rough 
surface  would  also  indicate  that  the  isotropic 
assumption  is  also  incorrect.  To  date  we  have  not 
quantified  the  impact  of  these  two  incorrect  assumptions 
upon  the  results  to  be  shown  shortly. 

The  perturbation  result  for  Model  III  (no  layer. 
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roughness  at  the  pressure  release  boundary) ,  is  well 
known  and  given  by  the  double  integral  on  the  first  line 
of  Fig.  5.  Theta  sub  zero  is  the  incident  angle,  small 
k  perpendicular  and  capital  k  perpendicular  are  the 
horizontal  wave  number  components  of  the  incident  plane 
wave  and  the  rough  surface  respectively.  P5  is  the  two 
dimensional  spectrum  of  the  rough  surface  which  can  be 
related  to  the  measured  one  dimensional  spectrum  *p[,  if 
we  assume  isotropic  roughness. 

We  fit  the  measured  data  for  Pj'  with  the 
analytical  expression  given  in  the  middle  of  the  third 
line  up  from  the  bottom  and  find  that  it  depends  on 
sigma,  the  correlation  length  and  goes  like  F  to  minus 
three  for  large  values  of  spatial  frequency.  Using  the 
analytical  fit  instead  of  the  actual  data  proves 
advantageous  for  we  can  then  find  an  analytic  expression 
for  the  one  dimensional,  auto  correlation  function,  MS) 
which  is  given  in  terms  of  the  correlation  length  and1  a 
modified  Bessel  function  of  the  first  kind. 

We  also  can  find  the  two  dimensional  isotropic 
spectrum,  which  has  a  form  similar  to  the  one 
dimensional  spectrum  except  that  it  behaves 
asymptotically  like  F  to  the  minus  four. 

Finally,  after  some  intricate  but  straightforward 
manipulations,  the  scattered  reflection  coefficient  is 
given  by  the  integral  on  the  bottom  line.  This  integral 
is  evaluated  numerically  for  various  frequencies  and 
angles  of  incidence. 

Two  approximations  (Fig. 6)  to  the  scattering 
integral  of  Fig.  5  have  been  f requently^sed  in  the 
past.  For  the  Kirchhoff  approximation  ^  is  replaced  by 
a  delta  function  which  yields  an  auto  correlation 
function  approximately  equal  to  unity  and  thus  ignores 
any  contributions  from  the  slopes  of  the  rough  surface. 
The  reflection  coefficient  is  then  found  to  have  a 
squared  dependency  with  both  frequency  and  the  cosine  of 
the  incident  angle. 

The  Marsh  approximation  which  would  appear  to  be 
valid  for  shallow  grazing  angles,  is  given  by  an 
abbreviate  version  of  the  scattering  integral. 

Utilizing  the  analytical  fit  to  the  one  dimensional 
spectrum  and  then  carrying  out  the  integration  one  finds 
that  the  reflection  coefficient  under  the  Marsh 
approximation  has  a  three  halves  dependency  bn  frequency 
a  linear  dependency  on  the  cosine  of  the  incident  angle 
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and  a  dependency  on  the  correlation  length. 

.  The  reflection  loss  for  all  three  results  namely 
the  exact  MSP  integral  shown  as  the  solid  line,  the 
Kirchhoff  approximation  given  by  the  line  with  long  and 
short  dashes  and  the  Marsh  approximation  given  by  the 
dashed  line  is  provided  in  Fig.  6  for  47  Hz  ando*  2 
and  1  =  44m.  For  incident  angles  between  77  and  90 
degrees,  which  is  the  region  of  interest  to  low 
frequency  Arctic  acoustics,  the  Marsh  approximation 
yields  a  greater  loss  and  the  Kirchhoff  approximation 
less  loss  than  the  exact  answer.  The  juxtaposition  of 
the  three  results  in  this  region  remains  the  same  for 
all  frequencies  below  100  Hz, 

Since  the  Marsh  approximation  yields  the  greatest 
reflection  loss  an  indication  of  how  well  the 
propagation  loss  predictions  from  all  three  results 
would  compare  with  the  empirical  result  can  be  obtained 
by  merely  examining  the  result  using  the  Marsh 
approx imation. 

The  FFP  prediction  of  propagation  loss  at  47  Hz 
using  the  Marsh  approximation  and  «  =  2m  and  correlation 
length  equal  to  44m  which  were  the  average  values  of  the 
ice  statistical  parameters  is  given  by  the  wavy  line  at 
the  top  of  Fig. 7.  The  empirical  fit  to  the  data  is 
provided  by  the  solid  line.  It  is  clear  that  the  Marsh 
approximation  provides  too  little  loss  at  this  frequency 
when  the  average  values  for  the  ice  statistical 
parameters  are  used.  This  discrepancy  persists  for  all 
frequencies  spanned  by  the  data. 

If  concurrent  measurements  were  not  available  and 
values  for  free  ice  statistical  parameters  could  be 
assigned  arbitrarily  perfect  agreement  with  the  data  can 
be  easily  attained  at  all  frequencies  with  the  Marsh 
approximation  merely  by  changing  G  to  3m  which  is  not  an 
unreasonable  value.*  The  result  at  47  Hz  (displaced  by 
20  dB)  is  shown  on  the  lower  portion  of  Fig.  7.  One 
could  also  obtain  excellent  agreement  with  the  Kirchhoff 
approximation  at  all  frequencies  if  <5  is  increased  to 
5m.  Thus  obtaining  agreement  with  experimental  data  in 
the  absence  of  concurrent  measurements  or  equivalently 
by  using  free  parameters  is  felt  to  be  a  poor  criteria 
for  judging  the  validity  of  scattering  formulas. 

If  the  scattering  surface  of  Model  III  were  taken 
as  the  ice-air  interface  instead  of  the  water-ice 
interface  we  suspect  the  comparison  using  the  average 


*  The  measured  values  of  O  varied  between  1  and  3m. 
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values  for  ff  and  1  would  be  worse  since  &  for  the 
topside  roughness  is  a  smaller  owing  to  melting  and  the 
effects  of  wind  erosion.  We  thus  conclude  that  Model 
III  is  inappropriate. 

The  perturbation  answer  for  Model  II  (liquid  layer 
with  no  roughness  at  the  top  and  no  absorption  loss 
within  the  layer)  is  given  in  terms  of  the  unperturbed 
reflection  coefficient  R^  and  the  double  integrals  I  sub 
one  and  I  sub  two  provided  in  Fig.  8.  The  M  title 
functions  are  involved  expressions  for  the  acoustic 
fields  in  each  region  evaluated  at  the  mean  surface  and 
are  provided  in  Fig.  9. 

We  make  two  observations  which  hopefully  provide 
some  insight  into  the  results  for  Model  II.  First  the 
unperturbed  reflection  coefficient  has  unit  magnitude 
for  all  incident  angles.  This  results  because  the 
absorption  in  the  layer  has  been  set  to  zero  and  the 
pressure  release  surface  behind  the  layer  perfectly 
reflects  any  transmitted  energy  back  into  the  water. 

Thus  the  presence  of  a  critical  angle  has  no  special 
significance  to  the  problem.  The  second  point  to  note 
is  that  the  difference  in  densities  between  the  two 
media  is  almost  zero. 

The  reflection  loss  for  Model  llf  under  the 
Kirchhoff  approximation  and  with  <7  =  2  and  1  =  44m  is 
given  by  the  dotted  line  in  Fig.  10.  Recalling  that  the 
unperturbed  reflection  coefficient  has  unit  magnitude, 
the  influence  of  the  roughness  at  the  bottom  of  the 
layer  is  determined  by  the  deviation  of  the  curve  from 
zero.  The  Kirchhoff  result  for  Model  III  (roughness  at 
pressure  release  surface)  discussed  previously  is 
overplotted  as  the  solid  line.  Since  the  solid  line 
result  provided  too  little  loss,  utilizing  the  average 
statistics,  (see  upper  portion  of  Fig.  7)  then  it  is 
clear  that  the  liquid  layer  will  fall  short  by  a  wide 
margin.  As  frequency  increases  the  deviation  of  Model 
II  from  zero  does  increase  slightly  but  in  all  cases  is 
considerably  less  than  the  Kirchhoff  result  for  Model 
III.  We  believe  the  factor  controlling  the  scattering 
for  Model  II  is  the  difference  in  density  on  either  side 
of  the  scattering  surface.  Since  this  difference  is 
almost  zero  the  scattering  loss  is  small  and  it  is 
unlikely  that  the  answer  would  change  substantially  if 
the  full  MSP  integrals  were  evaluated  instead  of  making 
the  Kirchhoff  approximation. 

We  thus  conclude  that  Model  II  although  apparently 
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more  realistic  than  Model  III  actually  provides  poorer 
results  and  thus  is  also  inappropiate. 

This  summarizes  our  progress  to  date.  In  the 
future  we  shall  examine  the  answers  provided  by  Model  I 
which  was  the  full  elastic  layer.  *  Considering  the 
number  of  free  parameters  which  exist  for  that  problem, 
however,  we  do  not  believe  a  convincing  case  can  be  made 
until  measured  values  for  the  elastic  properties  of  sea 
ice  become  available.  We  also  plan  to  examine  the 
influence  of  a  rough  surface  which  is  not  isotropic  and 
where  statistics  are  non-stationary. 


*  Subsequent  to  this  presentation  W.  Kuperman  of  NRL 
and  H.  Schmidt  of  SACLANT  RNN  Model  I  using  realistic 
values  for  shear  speed,  absorption  within  the  ice,  and 
the  measured  roughness  parameterms  of<x=2m  and-£  =4m. 
The  scattering  loss  using  the  full  elastic  layer  was 
found  to  be  comparable  to  that  obtained  from  Model  II 
and  substantially  less  than  that  from  Model  III  which 
was  substantially  less  than  the  measured  scattering 
loss . 
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MATHEMATICAL  MODELS  FOR  SCATTERED 
PLANE  WAVE  REFLECTION  COEFFICIENT 

MODEL  I  (FULL  ELASTIC  LAYER) 

Z,  =  0  (PRESSURE  RELEASE) 

Pi.  (C,)p.  (Ci)s.  (ttl)p,  (£*i)s  Ojt  III  Obi 

Z2=h 

/  \  P*>  ^2  >  ' 

MODEL  II  (LIQUID  LAYER) 

VI/YvVWWVtAV^  Z1  =  o  (PRESSURE  RELEASE) 

Pi,  Ci,  Off,  aj,  *t,  (Jq,  Ib 

Zj=h  ' 

r|  \  p2.  C2 

MODEL  III  (NO  LAYER,  ROUGHESS  AT 
PRESSURE  RELEASE  SURFACE) 

Pi  =  0  Gj  OR  (7g,  tj  OR  /g 

Zi=° 

R  \  P2,  C2 

•  START  WITH  MODEL  III,  WORK  UP  TO  MODEL  I 
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a  )  MODEL  III  (NO  LAYER,  ROUGHNESS  AT 
J/J  PRESSURE  RELEASE  SURFACE) 

f  FULL  PERTURBATION  THEORY 

2  ^  ^  1/2 

R(f,  % ,o,t)  =  -1+k  cos Bq  ~~~~  £ d  k* P2  (K* )  {k2  -  (k*-K*  )•  (k^-  K 

;  HORIZONTAL  WAVENUMBER  COMPONENT  OF  INCIDENT  PLANE  WAVE 
K*  ;  HORIZONTAL  WAVENUMBER  COMPONENT  OF  ROUGH  SURFACE 

A  I 

—  ASSUMING  ISOTROPIC  ROUGHNESS 

P^l  WH,CH  ,s  RELATED  TO  MEASURED  ? 

—  USING  ANALYTIC  FIT  TO  P^F),  F  SPATIAL  FREQUENCY  FIND 

~  7 r  (2a//)2 

N,  (p)  =  <2/0  K,  [(2//  )pj  P,(F)  =  [(2//)2  +  (2ir  p)2]  3/2  “  F  F°R  LARGE  F 

AND  2 

pi  (FI  - - 4*(2g/^  2  =  F*4  FOR  LARGE  F  * 

21  [<2/02  +  <2itF)21 

2 

w/2  cos2  0  sin  0  |l+(  -”)  (sin2  0+sin2  0q)| 

R=-1+k4  cos  dQ  (a  l  )2 J  — - - - 2 - - - 2 - ~2  3/2 

{  [1+ (  y)  <«"2  0+sin2  *o  ”  [(y)  $in  a  $in  ®  o]  } 
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MODEL  III  (NO  LAYER,  ROUGHNESS  AT 
PRESSURE  RELEASE  SURFACE) 

KIRCHHOFF  APPROXIMATION 


k(>  1.  M  CO!  e0  >  1  (APPROPRIATE  FOR  STEEP  GRAZING  ANGLES) 


«v  ■»  o  -*  dNo 

P2  (K_l>  as  2ir  (K^)  N2  (p)as  1  (p)  a  0 

R  (f  ,fl0,  o)  as  -1+2  a2  k2 cos2  »0 


MARSH  APPROXIMATION 

kl>  1,k/cos200  <1  (APPROPRIATE  FOR  SHALLOW  GRAZING  ANGLES) 

I 

R  ^-1+2^2  k3/2  cos  eQ  o2f  dlKjJPg  (IKJ)  |KA|3/2  , 

USING  MEASURED  FIND 

.3/2 

■v  ,  k  9  iir/i 

R  (f,  e0,o,t)2z  -1+2.396 cos0Q  o*  e,7r,H 
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Figure  5 


REFLECTION  LOSS  (dB 


MODEL  ill  REFLECTION  COEFFICIENT 
COMPARISON  FREQUENCY  47  Hz 
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PROPAGATION  LOSS  (dB) 


COMPARISON  OF  MARSH  APPROXIMATION 
WITH  EXPERIMENTAL  DATA 

FFP  WITH  MARSH  APPROXIMATION 

- FIT  TO  EMPIRICAL  DATA 

F  =  47  Hz  a  =  .04387  dBfkyd 


/ 
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Figure  7 


MODEL  II  (LIQUID  LAYER,  WITH  NO 
ATTENUATION,  NO  ROUGHNESS  AT  TOP) 
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Figure  8 


MODEL  II  (LIQUID  LAYER,  WITH  NO 
ATTENUATION,  NO  ROUGHNESS  AT  TOP) 
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MODEL  II  REFLECTION  COEFFICIENT 

FREQUENCY  47  Hz 


INCIDENT  ANGLE  (0O) 


Figure  10 
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